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Abstract: Higher spin Dirac operators on both the continuum sphere(5 2 ) and its fuzzy 
■ analog^!*) come paired with anticommuting chirality operators. A consequence of this is 

seen in the fermion-like spectrum of these operators which is especially true even for the 

f . v ,, 

^-j- . case of integer-spin Dirac operators. Motivated by this feature of the spectrum of a spin 

, 1 Dirac operator on Sp, we assume the spin 1 particles obey Fermi-Dirac statistics. This 



o 



0> choice is inspite of the lack of a well defined spin-statistics relation on a compact surface 

such as S 2 . The specific heats are computed in the cases of the spin ^ and spin 1 Dirac 
operators. Remarkably the specific heat for a system of spin | particles is more than that 
of the spin 1 case, though the number of degrees of freedom is more in the case of spin 
1 particles. The reason for this is inferred through a study of the spectrums of the Dirac 
operators in both the cases. The zero modes of the spin 1 Dirac operator is studied as a 
function of the cut-off angular momentum L and is found to follow a simple power law. 
This number is such that the number of states with positive energy for the spin 1 and spin 
^ system become comparable. Remarks are made about the spectrums of higher spin Dirac 
operators as well through a study of their zero-modes and the variation of their spectrum 
with degeneracy. The mean energy as a function of temperature is studied in both the 
spin i and spin 1 cases. They are found to deviate from the standard ideal gas law in 2+1 
dimensions. 
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1. Introduction 

The fuzzy 2-sphere(S'^) fl|, |2[ is an example of a noncommutative spacetime where the 
algebra of functions on the commutative 2-sphere(S' 2 ) is approximated by an associative 
matrix algebra. In principle this helps in regularizing field theories on such spaces, providing 
an alternative way to discretize spacetimes. The most non-trivial feature of this approach 
is the preservation of symmetries of the commutative spacetime even at the discrete level. 
Field theories on Sp have been actively pursued in the last few years || ||, [5], ||. Several 
works on numerical simulations of scalar fields and gauge fields on Sp have also been 
done. |8|, ^, [U], 11, |l!|. These works study phase transitions of scalar fields on the 



fuzzy sphere and they show the existence of new phases in the continuum which break 
translational as well as other global symmetries. The study of a real scalar field with a 
(p 4 interaction reveals the existence of a non-uniform ordered phase which is similar to the 
striped phase 13, 14, |D|. This feature distinguishes it from its commutative counterpart. 



Apart from these novel phenomena one can also incorporate supersymmetry in a precise 
manner on Sp [16]. 



Fuzzy spaces other than Sp have also been seen in the context of topology change [17 



|18|| . Physics on other fuzzy spaces have also been considered [19, po| [. 

To add to these new features on Sp, in this paper we probe the thermodynamics of 
spin 1 and spin ^ particles on the fuzzy sphere. We find several counterintuitive results 
which we will present in this work. We work with the spectrum of the Dirac operators for 
the spin 1 and spin \ case. This is natural to do as the Dirac operator is fundamental to 
physics and is useful in formulating metrical, differential geometric and bundle-theoretic 
ideas. Moreover in Connes' approach to noncommutative geometry [plf , the Dirac operator 
gains fundamental significance as part of the spectral triple in formulating the spectral 



action principle [22]. There are several ways to construct the Dirac operator on the fuzzy 
sphere |23|, | 



25, 26]. They construct the Dirac operator for a spin ^ particle. In our 



approach we construct the Dirac operators using the Ginsparg-Wilson(GW) algebra 
Through this method we can extend the construction to all higher spins as studied in [28]. 

We consider the spin 1 Dirac operators constructed in p3]. Unlike the spin | case, 
analytic computation of the spectrum of these operators are difficult. This lead us to 



compute its spectrum numerically [29]. Though we could not go to arbitrarily large values 
of the cut-off angular momenta L, we could still predict the spectrum's behavior in the 
continuum by observing the striking patterns that emerged for the spectrums for the values 



we could compute for. In 2J], 3 different Dirac operators along with 3 chirality operators 
were constructed for the spin 1 case. They are all unitarily inequivalent as proved in p9|[ . 
We consider only the spectrum of the traceless spin 1 Dirac operator in this work. 

Using this spectrum we first compute the partition function for a system of spin 1 
particles on Sp. For doing this we need to assume the particles obey a particular statistics. 
As we are dealing with a chiral system we assume that the particles obey the Fermi-Dirac 
statistics. However it should be noted that the conventional proofs of the spin-statistics 
theorem hold in relativistic quantum field theories (qft's) in three or more dimensions. 
They use the axioms of local relativistic qft's. For comprehensive proofs see [3l| ]. 
Field theory on the fuzzy sphere is not a relativistic one as the symmetry group of the 
underlying theory is SU(2). This being the case there is no well defined spin-statistics 
relation on the fuzzy sphere. However there are spin statistics relations which do not require 
relativity and which are topological [32, 33, 63]. General theory for quantum statistics 



in 2 spatial dimensions have also been discussed 133] . The non-triviality in two spatial 
dimensions arises due to topology of the configuration space of indistinguishable particles 
living on such a space. The fundamental group for such a configuration space( ^ R g A , 
where A = (r*i, f*2, fj^)\fi € M 2 and fj = fj for some i ^ j and Sn is the symmetric group 
of N particles.) is the braid group B^. For the case of S 2 instead of M 2 , the fundamental 
group is still the braid group with an additional constraint 36, 37]. These considerations 
allow for the possibility of the assumption of anyonic statistics |38|, 10] in our case, 
but we do not consider these possibilities in this work and only briefly remark about them 
in the final section. 

The main result of this paper is the mean energy of the spin ^ system is more than 
that of the spin 1 system. This is surprising given the fact that there are more number 
of states possible in the spin 1 case, 3 x (2L + l) 2 than in the spin i case, 2 x (2L + l) 2 . 
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Nevertheless we come to terms with this strange behavior by looking at the distribution of 
the eigenvalues of the Dirac operators in both the cases. In particular the spectrum of the 
spin 1 Dirac operator consists of a number of zero modes which are absent in the spin \ 
case. These characteristics provide an answer to the strange behavior of the mean energies. 
A consequence of this is also seen in the specific heats of the two systems. The entropy of 
the spin ^ system is also more than the spin 1 system. 

The other result is the deviation of the plot of the mean energy vs temperature from 
the corresponding curve for an ideal gas on a two dimensional space. The mean energy 
of an ideal gas of massless particles on a flat two dimensional space goes as T 3 . We find 
a deviation from this law which is attributed to the dispersion relation of the spin \ and 
spin 1 systems on the 2-sphere. 

The paper is organized as follows. Section 2 reviews the noncommutative algebra of 
functions on Sp. The GW algebra and the construction of the Dirac operators are reviewed 
in section 3. In section 4 we describe a way to find the spectrum of the spin 1 Dirac operator 
for arbitrarily large cut-off L. This is done by looking at a particular scaling behavior found 
in their spectrum. The grand canonical partition function is computed in section 5. This 
is used to compute the mean energy and the specific heats for both the cases of spin 1 and 
spin ^ systems. This section presents the relevant numerical results. Zero mode analysis 
is also carried out in this section along with the reasons for the strange behavior. We also 
speculate the behavior of the mean energies for higher spin Dirac operators. Section 6 
discusses the deviation from ideal gas behavior of these systems. We conclude in section 7 
with a few remarks. 



2. Geometry of Sp 

The algebra for the fuzzy sphere is characterized by a cut-off angular momentum L and is 
the full matrix algebra Mat(2L + 1) = M 2 l+i of (2L + 1) x (2L + 1) matrices. They can 
be generated by the (2L + l)-dimensional irreducible representation (IRR) of SU (2) with 
the standard angular momentum basis. The latter is represented by the angular momenta 
Lf acting on the left on Mat(2L + 1): If a £ Mat{2L + 1), 

Lfa = Ua (2.1) 

[Lf,Lf]=ie ijk Lf (2.2) 

{Lff = L(L + 1)1 (2.3) 

where Li are the standard angular momentum matrices for angular momentum L. 
We can also define right angular momenta Lf: 

Lfa = aLi,a£ M 2L +i (2.4) 

[Lf,Lf] = -ie ijk L* (2.5) 
(Lff = L{L + 1)1 (2.6) 



We also have 

[Lf,L?}=0. (2.7) 

The operator £j = Lf — Lf is the fuzzy version of orbital angular momentum. They 
satisfy the SU(2) angular momentum algebra 

[£i,Cj] = ieijkCk (2-8) 

In the continuum, S 2 can be described by the unit vector x 6 S 2 , where x.x = 1. Its 
analogue on Sp is ^j- or such that 

lim ^— = x,. (2.9) 

This shows that do not have continuum limits. But £j = Lf — Lf does and becomes 
the orbital angular momentum as L — > oo: 

lim Lf - Lf = A (2.10) 

3. Construction of the Dirac Operators 

In algebraic terms, the GW algebra A is the unital * algebra over C , generated by two 
^-invariant involutions r,r'. 

A = {r, r' : r 2 = r' 2 = i , r* = r , r'* = r'} (3.1) 

In any * -representation on a Hilbert space, * becomes the adjoint f. 
Consider the following two elements constructed out of T,r': 

r 1 = ±(r + r'), (3.2) 

Y 2 = l -{T-T'). (3.3) 

It follows from Eq.([0]) that -{Ti,^} = 0. This suggests that for suitable choices of T, 
r', one of these operators may serve as the Dirac operator and the other as the chirality 
operator provided they have the right continuum limits after suitable scaling. 

For the spin 1 case the combination which leads to the desired Dirac and chirality 
operators were found in [28] and they are 

D\ = L p^; 1 ^ 1 ) , (3.4) 
D 2 = L p-i~ r £+^ (3.5) 

and 

L r R 



Dz = L{ ^ n L . (3.6) 



with 



71 = 2 ' (3 - 7) 



72 = ^ ^ J (3.8) 

and 

73 = (3-9) 
as their corresponding chirality operators. In the above equations 

r 2(t.L L + L + 1)(£.L L + 1) — (L + 1)(2L + 1) 



(3.10) 



i+1 (L + 1)(2L + 1) 

_ 2(-E.L R + L + l)(-g.L fi + 1) - (L + 1)(2L + 1) 
iL+1 " (L + 1)(2L + 1) ' ( - Kl " 



^ 2(£.L i -L)(£.L L + l)-L(2L + l) 



^ = — L(2L + 1) ^ 



B 2(E.L fi + L)(E.£ fl - 1) - L(2L + 1) 
L(2L + 1) 



(3.13) 



L _ -2(EX L - L)(E.L L + L + 1) - L(L + 1) 
Fi " L(L + 1) ' (3 - M) 



and 



(3.15) 



R = 2(Z.L R + L)(-E.L R + £ + 1) - L(L + 1) 

L " L(L + 1) 

The operators in Eq.( [3.10 )-Eq,( |3TT5| ) are generators of GW algebras and are obtained from 
left and right projectors to eigenspaces of the total angular momentum, L + £, where E 
are the matrices representing the spin 1 representation of SU(2). 
The continuum limits of Eq.Q3.4|)-Eq,(|3.6|) are 



Dx = (£X - (£.x) 2 + 2) + 2(£.x) + {EX, £.x}, (3.16) 

D 2 = (EX - (£.x ) 2 + 2) - 2(£.x) - {EX, S.x } (3.17) 

and 

X) 3 = EX - (S.x) 2 + 2. (3.18) 
The corresponding chirality operators in the continuum are 

71 = (£.x) 2 + (£.x) - 1, (3.19) 

72 = (E.x) 2 - (£.x) - 1 (3.20) 

and 

73 = 1 - 2(S.x) 2 (3.21) 

respectively. 



4. Spectrum of the Spin 1 Dirac operator for large N 



For a given cut-off L = (N — l)/2, the eigenvalues of the spin 5 Dirac operator is given by 
±(j + I) where j is the eigenvalue of the total angular momentum. The degeneracy (g) of 
each of these eigenvalues (E) is given by 

g = 2j + 1. (4.1) 

This gives the relation between the energy eigenvalue and its degeneracy as 

Ei = ±-. (4.2) 

In contrast, the spin 1 system is richer with varied eigenvalues for each L including a 
number of zero modes. Though we cant find the corresponding relation between the energy 
eigenvalue E and its degeneracy g analytically we can do so numerically | |29f . 

The numerical results show that g increases with E for small values of E. After 
reaching a maximum for an intermediate value of E, g decreases with E. This is unlike 
the spin ^ case where g is linear in E. After rotation (around origin, by a fixed angle) the 
plot of E' vs g' fits perfectly with a parabola. The rotation angle was found to be cut-off 
independent. Encouraged by this we tried to find if there is any universality in E' vs g'. 
We found that results for different cut-off lie on a universal curve. This is done by scaling 
E' and g' so that the minimum of different curves match. In figure. 1 one can see the results 
for N = 21,27 and 45 lie on top of each other after scaling. This result strongly suggest 
that we can obtain the spectrum for arbitrarily larger cut-off (N) if we know the scaling 



variables for different N [29]. The universality also suggests that E' vs g' is described by 

only two independent parameters. Once we know E' vs g' we rotate it to get the relation 
between E and g, which is given by 

y/b( g y-4ac(g)-b(g) 

El = 2a (43) 

where 

a = acos 2 9, (4.4) 

b{g) = 2a cos 6 (sin Og + 77) + sin 6, (4.5) 

c(g) = a(sin 6g + r/) 2 + (3 — cos 6g (4.6) 

and 

9 = 2.26159 radians. (4.7) 

As discussed above only f3(N) and r](N) are independent parameters. For a we need 
to know its value for a particular N. For all other values of N, a can be calculated using 
P(N) and rj(N). 
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Figure 1: Scaling property of the spectrum of the spin 1 Dirac operator. 
5. Numerical Results: I 

5.1 The Partition Function and the Mean Energy 

As explained in the introduction, we assume the spin 1 particles to obey fermionic statistics. 
The grand-canonical partition function is given by 

InZ = ^ 9%ln (l + e'^-^A (5.1) 
-i=i 

where gi is the degeneracy of the ith level, E% is the energy of the ith level, \i is the chemical 
potential and /? = ^— ^ . In the commutative case when there is no cut-off, the product over 
i extends till infinity, but here we are restricted by the cut-off angular momentum L. 

For the spin 1 case we numerically computed the spectrum of the Dirac operator 
given by Eq.(|3.6|) in |29[|. We do not know how to find an analytic expression for the 



spectrum and so we compute the partition function numerically. (See however [29| for an 
analytic expression for the spectrum of the spin 1 Dirac operator derived as a result of the 
numerical computations.) The analogous situation for the spin ^ case is fa better as we 
know its complete spectrum analytically for arbitrarily large cut-off L. 

From the grand-canonical partition function in Eq.([|[l|) we can use the standard for- 
mula to compute the mean energy which is 

dlnZ ^ gjEj 

{El) = -^= 2. e m-,) + i (5 - 2) 

i=l 

In what follows we take the Boltzmann constant fee = 1 and the chemical potential \x = 0. 
The above expression for the mean energy is used for both the spin 1 and the spin ^ cases. 
For the spin ^ case it becomes 

(El)= 2 yl + . (5. 3 ) 
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Figure 2: The mean energies of the spin 1 and spin \ systems. 



Note that in the above Eqns.(^2|, ^) the sums are restricted by the cut-off L. The mean 
energies for both the cases were computed for various temperatures from 0.1 to 50. We 
show only these plots here though we did go to higher values of temperature and found 
nothing new. The plot for the mean energies of both the spin 1 and spin ^ systems is 
shown in figure 2. The value of cut-off is L = 

In figure 2, the green curve shows the mean energy for the spin ^ system as a function 
of temperature and the red one shows the corresponding curve for the spin 1 system. The 
curves become flat for higher values of temperature. This is due to the presence of the cut- 
off angular momentum in our sum. If we go to higher values of temperature this flattening 
occurs towards the higher temperatures considered. The plot clearly shows that the mean 
energy of the spin \ system is much higher than the spin 1 system. This is inspite of the 
spin 1 system having more number of degrees of freedom than the spin | system. We know 
of no such analogous behavior in higher dimensions. 

Another interesting feature in the behavior of these curves is the crossing of the two 
curves for low values of temperature. This is not clear in figure 2 but is shown explicitly 
in figure 3. This plot shows that the mean energy of the spin | system is smaller than the 
spin 1 system till about T = 10.77 after which it stays above the spin 1 curve. 

We now try to explain the cause of this unusual behavior by looking closely at the 
distributions of the eigenvalues of the Dirac operators of the two systems. 
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5.2 Reasons for the strange behavior 

The main reason can be understood once we look at the spectrum of the Dirac operator in 
the two cases. 

Using the expressions for the energy as a function of the degeneracy we can study the 
differences between the two systems. The relation is given by Eq.(4.2) for the spin ^ case 
and Eq.( f4.3| ) for the spin 1 case. 
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Figure 3: The crossing of the mean energy curves of the spin 1 and spin § systems. 
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Figure 4: Energy as a function of the degeneracy g for the spin \ and spin 1 systems. 



This plot of the energy eigenvalues as a function of their degeneracies is shown in figure 



4. 



For a given cut-off L, figure 4 clearly indicates that eigenvalues of the spin ^ Dirac 
operator exceeds that of the spin 1 Dirac operator except for small values of the degeneracy 
g. The plot in figure 4 is shown only for positive values of the energy eigenvalue E. In 
the spin ^ case the energy eigenvalues linearly increase with the degeneracy g and so the 
maximum eigenvalue occurs for the j value 2L — 5 . We have ignored the maximum j value 
of 2L + i as they correspond to unpaired eigenstates of the Dirac operator, which will be 
inconsistent given the chirality of the spin | system squares to 1. This is a feature of the 
operator a.C in the spin ^ Dirac operator which is 



Di_=aX + -, 
2 2 



(5.4) 



where a are the Pauli matrices. 

In the spin 1 case the maximum eigenvalue is L and this occurs for some intermediate 
value of the degeneracy g as can be seen in figure 4. The reason why the spin 1 Dirac 
operator consists of all eigenvalues ranging from —L to +L for a given L can be seen by 



looking at the operator in the continuum given by Eq.( 3.18| ) which we write here again: 

D 3 = EX- (S.x) 2 + 2 (5.5) 

where £ are the matrices of the spin 1 representation of SU(2). The term S.x makes the 
analytic computation of the spectrum in the spin 1 case difficult when compared to the 
spin g case. We believe this term to be also the cause of the varied spectrum of the spin 1 
system. 

5.3 Zero modes of the spin 1 Dirac operator 

The spectrum of the spin 1 Dirac operator consists of a number of zero modes for each 
cut-off angular momentum L. The number of such zero eigenvalues follows a simple power 
law as a function of L. This number was found to be (2L + l) 2 + 2. This is an exact result 
and can be found analytically as explained in |2£| . This has also been verified numerically. 
With this result it follows immediately that the number of positive eigenvalues of the spin 
1 system are (2L + l) 2 — 1. The spin ^ Dirac operator has no zero modes as it has non-zero 
trace. Removing the states corresponding to the top mode gives us 4L 2 + 2L states with 
positive eigenvalues. The zero modes in the spin 1 case drastically reduce the total number 
of states corresponding to positive eigenvalues to (2L + l) 2 — 1 but this is still more than 
the corresponding number of states in the spin ^ case. 

The counting of the zero modes and the behavior of the spectrum with degeneracy in 
the two cases justify the counter-intuitive behavior of the mean energies. 

We now digress a bit to remark about the plot in figure 4. We try to speculate the 
energy versus degeneracies curves for higher spin Dirac operators. To do this we first find 
the number of zero modes for higher spin Dirac operators. We will compute this for the 
integer spin case. 

To construct higher spin Dirac operators on S F , we need to construct operators acting 
on Mat(2L + 1) ® C 2fc+1 where k is the desired spin. The spectrum of these operators will 
in general be hard to compute due to the presence of £.£ terms just as in the spin 1 case. 

The analytic computation of the number of zero modes was given in |29|. We extend 
those arguments to higher spins in the following. Consider the spectrum of the total angular 
momentum j for a given cut-off L: 

Spec J G {0, 1, 2, • • • , 2L - k, ■ ■ • , 2L, • • • ,2L + k}. (5.6) 

For an even-integer spin k, the number of zero-modes can be found by computing the 
following sum 

2L-fe | 

£(2i + l)+ [2(2L + 2j) + 1] = (2L + l) 2 + k 2 + k. (5.7) 



For an odd-integer spin fc, this number is 



k — 1 

2L-k-\ — 



Y (2j + 1) + t 2 ( 2L + 2 -? + !) + !] = ( 2L + x ) 2 + fc2 + k - ( 5 - 8 ) 



7=0 ■•— + 1 

J - 2 



These computations hold as there exists a traceless Dirac operator for all integer spin Dirac 
operators on Sp. This is because we can construct an integer spin Dirac operator from the 



following combination of generators of GW algebra [28|: 



D k = L (p^j • (5.9) 

As tr(T L ) = tr(r^), this operator is traceless. 

In the case of the Dirac operators for half-integral spins, there exists no such combina- 
tions of generators of GW algebras which have trace. This makes the number of states 
with positive energy eigenvalues for spin k and spin k — i, for integer k, comparable. 

We can then go on to compute their mean energies and compare them. We suspect 
{E k _i) > (Ek) to hold but we have no analytic proof for this. We could however compute 
the spectrums for the two Dirac operators numerically and carry out this comparison, but 
we do not do this here and save it for future work. 

The reason why this is interesting is the following. It seems from the plot in figure 4 
that the behavior of for small values of g is similar for higher spins as well. We leave 
this as a conjecture as we have no analytic proof for this but do have strong reasons to 
suspect so. 

It is also very likely that the plots of E versus g for higher spins will fall below the 
E = | curve. This is expected due to the fact that higher spin Dirac operators contain S.x 
terms along with £.£ [p8| . £ is the 2k + 1 dimensional representation of SU(2) for some 
spin k. The T*.x term disrupts the linearity between the energy and degeneracy. It is easy 
to see this as a linear relation between the energy and the degeneracy is only possible for 
a Dirac operator which has just a £.£ term apart from constant terms. This can be seen 
analytically for any given spin k by looking at the spectrum of £.£: 

Spec £.£ = (k - m)(2j - k + m) - k 2 - m m E {0, 1, • • • , 2k}. (5.10) 

Only the spin i Dirac operator contains just the a.C term leading to the linear relation 
between its energy and their multiplicities. 

The Ti.x terms are present in the spin 1 case and it was remarked that these terms 
cause the energy versus degeneracy curve in figure 4. As these terms also occur for higher 
spin Dirac operators we expect a similar behavior from these systems. The reason why 
they occur for all higher spin Dirac operators is because of the fact that the fuzzy versions 
of these higher spin Dirac operators contain terms of the form 

{t,.L L ) n - (t.L R ) n 
L n_1 ' 




p 2 

Figure 5: Ratio of specific heats 

The continuum limit of these terms contain S.x terms. This is explained in detail in 

The preceding statements prove the non-linearity between the energy and their degen- 
eracies for all Dirac operators other than the spin ^ system. They however do not show 
that these curves fall below the corresponding curve for the spin ^ system. A complete 
answer to this question would only come from a numerical analysis of this system and at 
present we leave this question as a worthy one to explore in the future. 

5.4 Specific heats of the two systems 

The specific heat is defined as the derivative of the mean energy with respect to tempera- 
ture. A straightforward computation gives the specific heat as 



1 giEfeP Ei 



i ' 

As expected here too we find the specific heat of the spin \ system to be more than that 
of the spin 1 system. This is shown in figure 5. 

There is a region till T = 10.77 where the specific heat of the spin 1 system is more 
than that of the spin ^ system. This can be seen as a result of the crossing of the mean 
energy curves for the two systems as shown in figure 3. 

5.5 Entropies of the two systems 

The entropy is given by the equation 



This follows from 



S = InZ + P(E). 



(5.13) 



From these formulas it can be easily seen that the entropy of a spin 1 system is less than 
that of a spin \ system. 

6. Numerical Results: II 

6.1 Deviations from the Ideal Gas Law 

For a system of non-interacting massless particles obeying Fermionic statistics, the mean 
energy goes as T in 3+1 dimensions. This can be seen as follows: 



This makes the sum 



4L+2 
^2 T 9 



{E k ) = T £ TTT- (6 - 5) 

* 1 e 2 + 1 



x= 



As the limits of the sum depend on the temperature T we get no definite relation between 
the mean energy and temperature. It should be noted that the upper limit is dependent 
on T due to the cut-off L. We can remove this by allowing L to go to oo. In such a case, 
as already mentioned the sum becomes an integral making the dependence go as T 3 . This 
still does not remove the T dependence from the lower limit of the integral. This is due to 
the dispersion relation for the spin ^ particle which goes as j + \ . The additional \ can be 
attributed to the curvature of the sphere the system lives on. This results in the deviation 
from ideal gas law on 2+1 dimensional space. 



DO 
O 



where p is the energy of the massless particle. This is the dispersion law for a massless 
particle on a flat space which has the Poincare group has its group of symmetries. We now 
substitute 

V 

X= T 

to find 

/•oo 3 

(E)<x T*j o cfc— . (6.2) 

In a similar manner, in 2+1 dimensions, the mean energy goes as T 3 . This law holds 
however only for a system living on a flat 2+1 dimensional spacetime. 

In the case of the spin \ system living on S 2 or Sp the mean energy is given by 

2L+- 

J^o 2 e^T~ +1 

We have cut-off the sum with a cut-off L. If we arbitrarily increase the value of the cut-off 

L we will find the sum replaced by an integral over j and the upper limit in the sum goes o 

to oo. In the above equation make the substitution 

^±1 = x. (6.4) 
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Figure 6: Deviation from ideal gas law for the spin \ system. 
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Figure 7: Deviation from the ideal gas law for the spin 1 system. 



Similar arguments hold for the spin 1 case also. This can be easily seen from our 
analytic expressions for the spectrum of the spin 1 Dirac operator as seen in the previous 
section. We do not write the simple details of this here. 

The deviations for the spin ^ and the spin 1 system are shown in the plots in figures 
6 and 7. 



7. Conclusions 

Spin systems on the fuzzy 2-sphere show several interesting features as discussed in this 
paper. The main result being the mean energy of a spin | system is greater than that of 
a spin 1 system. This is a novel phenomenon with no known higher dimensional analog. 
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Its implications are worth exploring. The simple power law of the zero modes is another 
remarkable feature of these systems. 

Speculations were made about the spectrums of higher spin Dirac operators on Sp. 
Though we could not arrive at a concrete result, we conjectured that it is likely that the 
mean energy of the spin i system will be more than that of the mean energies of other 
spins greater than ^. This gives an upper bound on the mean energy of these spin systems 
on Sp. This bound will also translate to the entropies of these systems, thus making the 
entropies of these spin systems to be bounded by the entropy of the spin | system. This 
would mean that an increase in the number of states on Sp have no effect on the entropy, a 
property very much reminiscent of the holographic principle though there is no bulk theory 
involved here. This is conjectural and we plan to explore this in the future. 

The deviation from the ideal gas law is another interesting phenomenon in this fuzzy 
system. It can definitely serve as a model for explaining 2 dimensional systems which show 
such behavior. 

We carried out similar computations by assuming other exotic statistics like anyonic 
statistics, but found nothing interesting. For this we used the following formula instead of 
Eq.Q 



where a is allowed to vary from — 1 . . . , 0, • • - + 1 going through all possible statistics between 
bosons and fermions. However there are different approaches in formulating the statistical 
mechanics of anyons and our treatment is by no means complete. We do plan to study this 
further in the future as we maybe missing some connection. (See also |^2| in this regard.) 

Dirac operators appear in 2-dimensional graphene systems in condensed matter physics. 
It would be interesting to know if the continuum limits of our models will be useful for 
similar carbon systems which have a spherical shape like fullerene (Cqq). We do not know 
of any direct application now. 

Dirac operators can be constructed on other fuzzy spaces as well. It is a very interesting 
problem to see if the method of GW algebras can be used to construct higher spin Dirac 
operators on other fuzzy spaces. A study of their thermodynamics could lead to a lot of 
rich and strange phenomena like the ones found in this paper. 

These can have an application to cosmology in terms of the thermal history of the 



universe. We will report these results in a forthcoming paper [41|. Also see in this regard the 



paper [43] where an interesting toy model of the universe as a fuzzy sphere was presented. 

The new results found here are encouraging, and a hunt for applications of these 
features of fuzzy systems is a worthy endeavor. 
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